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ABSTRACT 

The mysterious 21 fim emission feature seen in sixteen C-rich proto-planetary nebulae 
(PPNe) remains unidentified since its discovery in 1989. Over a dozen of materials are 
suggested as the carrier candidates. In this work we quantitatively investigate eight 
inorganic and one organic carrier candidates in terms of elemental abundance con- 
straints, while previous studies mostly focus on their spectral profiles (which could be 
largely affected by grain size, shape and clustering effects). It is found that: (1) five 
candidates (TiC nanoclusters, fullerenes coordinated with Ti atoms, SiS2, doped-SiC, 
and Si02-coated SiC dust) violate the abundance constraints (i.e. they require too 
much Ti, S or Si to account for the emission power of the 21 /mi band, (2) three candi- 
dates (carbon and silicon mixtures, Fe2C>3, and Fe 3 04), while satisfying the abundance 
constraints, exhibit secondary features which are not detected in the 21 fj,va sources, 
and (3) nano FeO, neither exceeding the abundance budget nor producing undetected 
secondary feature s, seems to be a viable candidate, supporting the suggestions of 
iPosch et al.l (|2004f ). 

Key words: (stars:) circumstellar matter - infrared: stars - stars: AGB and post- 
AGB - stars: individual (HD 56126) 



1 INTRODUCTION 



Since its first detection in 1989 l|Kwok et al.lll989l b the so- 
called "21 fim feature" has been identified in sixteen proto- 
plane tary nebulae (PPNe; iKwok et al. l ll999l . lHrivnak et al.l 
2009) [and arguably also in two planeta ry nebulae (PNe ) 
associated with Wolf-Rayet central stars IIHonv et al.l 2001) 
and in two highly evolved carbon stars (|Volk et al.l 120001 )]. 
This feature has little shape variation among different 
sources with a peak wavelength at ~ 20.1 fim and a FWHM 
of ~ 2.2-2.3 fim. Most of these sources exhibit quite uniform 
characteristics: they are metal-poor, carbon-rich F and G 
supergiants with strong i nfrared (IR) exces s and over abun- 
dant s-process elements (IZhang et al.ll2006h . 

After its discovery, over a dozen of carrier can- 
didates have been proposed (see Fig. [TJ, with the 
number of proposed candidate materials comparable to 
the total number of the 21 fim. fe ature sources (see 
I Andersen et al.ll2005l . lzhang et all2006h . These include both 
inorg anic materials: (a) Ti C nanocl u sters l|Von Helden et all 
1200 Jl, (b) SiS? grains dGoebell Il993l b (c) Doped-SiC 



dust (|Speck fc Hofmeisterll2004r ). (d) carbon and silicon mix- 
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tures (Kimura ct al. 2005) , ( e) SiC core-SiQ2 m antle grains 
|Posch et a!ll2004r), (fi Fe O l|Posch et al.ll2004r i, (g) Fe 2 3 , 
and (h) Fe 3 o71co3[l990); and organic materials: (i) large- 
cage carbon particle s (fullerenes) coordinat ed with Ti atoms 
(iKimura et al.l l2005b (j) urea or thiourea (|Sourisseau et all 
Il992h . (k) polycylic aromatic hydrocarbo n (PAH), and (h) 
hydrogenated amorph ous carbon (HAC) l|Buss et al .11 19931 : 
lJusttanont et al.lll996lb 

How ever, the carrier o f this feature remains unidenti- 
fied (see lPosch et al.ll2004l for an excellent overview). It is 
considered as one of th e most interesting unresolved myster- 
ies in astrochemistry l|Kwok et al.l 12002 ). Previous studies 
mostly rely on a comparison of the spectral profile of a can- 
didate material with the observed 21 fim emission feature. 
In this respect, astronomers often consider spherical dust 
using Mie theory. However, the spectral profile is expected 
to vary with grain shape, size, and the presence of voids (e.g. 
sec jHuffmarJll989l :lLr et al . 20081: IVoshchinnikov fc Hennind 
l200sh . Although models based on spherical dust may not fit 
the observed 21 fim emission feature, one can not rule out 
the possibility that a reasonably good fit may be obtained by 
fine-tuning dust shape and size distributions. Indeed, some 
studies have demonstrated that the spectral fit is improved 
with a continuou s distribution of ellipsoids (CDE; see e.g. 
IPosch et ai]|2004 ) for the dust shape. We should note that 
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even the failure of CDE models can not readily rule out the 
considered candidate material since the spectral profile could 
be further affected by many other factors (e.g. , clustering of 
individual grains, see IRouleau fc Martin|[l99ll ). 

In this work we attempt to constrain the nature of the 
21 /jm feature carrier by performing a systematic study of 
the validity of the above-listed candidate carriers. In view of 
the fact that a spectral profile "mismatch" could probably 
be improved by considering dust shape and size distribu- 
tions, clustering of individual grains, and the presence of 
voids, unlike previous studies, we do not completely rely on 
whether the interested candidate material produces a band 
profile closely matching the observed one. Instead, we take 
an alternative approach which at least complements previ- 
ous approaches (which are mainly based on spectral profile 
matching): we focus on the elemental abundance required 
to account for the total power emitted in the 21 /im feature. 
This approach is less affected by the (unknown) dust clus- 
tering and shape distribution and is therefore more robust. 

This paper is organized as follows: we first present in ij2] 
the general constraints on the validity of a proposed carrier. 
We then apply these constraints to the above-listed candi- 
date materials (S|. For this purpose, we choose the pro- 
totypical 21 fim feature source HD 56126 as a comparison 
basis. The main results are summarized in ijS] 



2 GENERAL CONSTRAINTS: BAND 

STRENGTH AND ABUNDANCE BUDGET 

For a candidate carrier to be a viable explanation of the 
21 /im feature, in addition to the close match with the 21 fim 
feature observed in PPNe, it must satisfy the abundance 
constraint, i.e., the candidate material must be abundant 
enough to account for the total power emitted in the 21 fim 
feature. Further, it should of course not produce any sec- 
ondary features which are not observed in the 21 /im sources. 

The 21 fim feature is one of the strongest IR dust fea- 
tures in C-rich evolved objects. The strongest 21 fim source 
(HD 56126) emits ~ 8% of its total IR power in this feature, 
while the 21 /im feature in other sources amoun ts to ~ 1— 
7% of their total IR power (|Hrivnak et alj|2000l ). That the 
21 /jm feature emits such a large energy would put stringent 
constraint on the abundance of the carrier, particularly, on 
the abundance of the relatively rare element in the suggested 
carrier material (e.g. Ti in TiC). 

Let E to t be the total power emitted in the 21 /im band, 
which is an observational parameter. It is related to the total 
mass of the carrier dust Md through 



M d 



Kabs(A) x AnB\{Td)d\ , (1) 



21 /^m band 



where the integration is over the 21 /im band (but with the 
continuum underneath the 21 /im feature subtracted), K a bs 
is the mass absorption coefficient (also known as "opacity") 
of the dust, i?A(2d) is the Planck function of a black body 
with temperature Td at wavelength A, and Td is the dust 
temperature. 

If the carrier material contains element X, the abun- 
dance of element X (relative to H) in the dust can be written 



[X/H] = 



nx M d //j d m H 



(2) 



M H /m H 

where Mh is the total hydrogen mass of the circumstellar 
envelope of the object, iiih is the mass of a hydrogen atom, 
/id is the molecular weight of the dust grain, and nx is the 
number of atoms per molecule for element X. With a good 
knowledge of the dust temperature and the mass absorption 
coefficient in the 21 um wavelength range of the carrier, one 
can therefore estimate the abundance of a typical element 
locked up in dust required to emit the observed power E to t 
in the 21 um band 

nx Etot 



[X/H] = 



W M u J 2lMmband ^bs(A) x 4nB x (T d )d\ 



(3) 



It can be seen that for a given E to t, the required abun- 
dance [X/H] is inversely proportional to the integral of ab- 
sorption coefficient K a bs and the dust radiation intensity 
B\(Td) over the 21 um band range. 

We approximate the mass absorption coefficient profile 
with a Drude function 



^abs(A) — 



4bs 27A 



(A-A2/A) 2 +7 2 ' 



(4) 



where A c 



20.1 um and 7 ~ 2.2 /im are respectively the 

nt _ 



peak wavelength and FWHM of the profile □ and Kjjj 
/21pm band K abs(A)dA is the mass absorption coefficient inte- 
grated over the 21 um band. The adoption of a Drude func- 
tion and the values of A p and 7 which match the observed 
feature profile mostly is the best condition for the absorption 
coefficient. The abundance under such condition is the min- 
imum. By substituting Eq. Q into Eq. © , the lower limit 
of [X/H] becomes 

Tlx Etot 



[X/H] 



8MdMH7«kbB/ 2 i, 



\(T d ) 



-dX 



(5) 



1 Mm band ( x _ A 2 j a) T + 7 2 

Note that the integration is continuum-subtracted. Since the 
Planck function B X (T) is a monotonically increasing func- 
tion of Td, the minimum abundance of element X, [X/H] min , 
becomes a decreasing function of Td and KabV For a fixed 
Etot, given a reasonable dust temperature range (usually 
~ 100-200K for the circumstellar envelope of post-AGB 



1 The Drude profile, closely resembling a Lorentz profile and 
having more extended wings than a Gaussian profile, is ex- 
pected for classical damped harmonic oscillators jLill200sh ■ With 
A p = 20.1/mi and 7 = 2.2/im, the Drude profile fits well the ob- 
served 21 fim emission feature except the blue- wing (which is a bit 
more extended than that observed), while the Gaussian profile is 
a bit too narrow in the red-wing. To be more physical, we should 
fit the observed 21 fim emission feature with K a bsM B\(Tj), 
a sum of products of a Drude mass absorption profile and Planck 
functions of a range of temperatures (since the carrier of the 
21 fim feature is expected to have a range of thermal equilib- 
rium temperatures in the 21 /mi-feature-emitting shell). Indeed, 
we find that with A p = 20.1/im, 7 = 1.85/im, and T fb 90 K, 
the product of a Drude profile and a single-temperature black- 
body KabsM B\{T) closely fits the observed 21 fim emission fea- 
ture (see the inserted panel in Fig. lh). In this work, we take 
A p = 20.1 fim and 7 = 2.2 fim in order to maximize the energy out- 
put in the 21 fim feature so as to minimize the dust abundance 
requirement. 
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stars) and k^ s from laboratory data, we can estimate the 
range of [X/H] . . For a viable candidate, the required mini- 
mum abundance [X/H] min must not exceed what is available. 

If K a bs(A) is measured in the laboratory, K a b s can be 
obtained from the direct integration of the measured profile 
of K a bs(A). In a few cases, Ac a bs(A) is not directly measured, 
we obtain n'^ s from the absorption cross section C a bs or the 
absorption efficiency Q a bs calculated through Mie theory 



/, • ., = / K abs (A)dA 

'21 fim band 

C a bs (l, A) 
M 



d\ 

\ ! 

21 fim band 

3 f Qabs (a, A) 

J 21 band u 



dX 



(6) 



where M and a are respectively the mass and radius of the 
dust grain (we assume that the dust is spherical), and p d is 
the mass density of the dust. 



count for the 21 /im feature 



(M d /M H ) r 



V* [ X / H ] mil >x 



i J21 im 



B X (T A ) 



-dA 



(8) 



Apparently, for a valid carrier candidate, the minimum re- 
quired mass {Md/M-g) mln should not exceed the maximum 
available mass (M d /MH) max . 

In Fig.HJwe plot (M d /M H ) max against 

K abs f° r t ne P r °" 

posed carrier candidates. We also plot (M d /MH) min as a 
function of k^ s for a range of temperatures T d . Note that 
the temperature at which (M d /M H ) max » (M d /M H ) min is 
the lowest temperature the dust should have. For exam- 
ple, at T w 880 K TiC nanoclusters have (M d /M H ) max « 
(M d /M H ) min . In order for (M d /M H ) min not exceeding 
(M d /M H ) max , TiC dust must have T d > 880 K, otherwise 
there is simply not enough dust material. 



3 THE TESTER: HD 56126 

A successful candidate carrier should be able to explain the 
observed 21 /im feature in all sources. A failure in a sin- 
gle source would be sufficient to rule out the candidate. To 
examine whether the carriers can account for the observed 
feature strength, we choose HD 56126, a prototypical 21 fim 
feature source, as the tester. 

HD 56126 (IRAS 07134+1005) is selected for the follow- 
ing reasons: (1) It is the strongest 21 fim feature emitter and 
its emission flux was accurately measured; (2) It is one of 
the best studied 21 fim feature s ources and it s basic param- 
eters are accurately determined. iHonv et alj (|2003h built a 
detailed dust radiative transfer model for the circumstel- 
lar envelope and deriv ed the composition and mass o f the 
dust shell of HD 56126. IVan Winckel fc Revniersl (|200Ch per- 
formed a homogeneous photospheric abundance analysis of 
HD 56126 (as well as several other 21 (im feature sources). 
We summarize in Table Q] the key relevant stellar and cir- 
cumstellar parameters which will be used in later analysis. 
The parameter Mh is not very certain in the range ~0.16- 
0.44 M Q . We generously take the high end (M H « 0.44 M©) 
for our analysis. Since the required minimum abundance 
[X/H] . is inversely proportional to Mh (see eq.|SJ, if a 
model based on the high end Mh ~ 0.44 M© already ex- 
ceeds the elemental budget it should certainly be rejected. 

By assuming that the X atoms available in the cir- 
cumstellar envelopes around the 21 fjxci feature sources (i.e. 
[X/H] ) are all depleted in the dust species proposed as a 
carrier candidate, we obtain the maximum dust mass (rela- 
tive to H) of the dust species containing the key element X 
from Eq.rj 

(Af d /M H ) max = fi d [X/HL/nx . (7) 

This is the upper limit on the amount of dust containing 
element X available to account for the 21 /im feature. For a 
given dust temperature T d , from Eq.[5]we obtain the lower 
limit on the amount of dust (relative to H) required to ac- 



4 ASSESSING INDIVIDUAL CARRIER 
CANDIDATES 

We assess proposed individual carrier candidates by exam- 
ining (1) whether they are capable of emitting the observed 
large amount of energy in the 21 ycm band without requiring 
more dust material than available, and (2) whether the can- 
didate carrier produces (undetected) secondary features. In 
view of the crucial role of elemental abundances in this as- 
sessment, we divide the carrier candidates into four groups: 
titanium-bearing, sulfur-bearing, silicon-bearing and iron- 
bearing grains. 

The 3.3, 7.7 and 11.3 /im PAH features and the much 
more promin ent 30 /im feature which i s generally attrib uted 
to MgS dust |Goebel fc Moseley|[l985l : IHonv et alj|2002h are 
seen in fourteen of the sixteen known 21 /mi feature sources 
l|Kwok et al.lll999l . iHrivnak et al.ll2009h . Except these fea- 
tures, the 21 /xm sources do not universally show any addi- 
tional dust features. However, in addition to the 21 /mi fea- 
ture, some of the proposed carrier candidates exhibit strong 
spectral features at other wavelengths as well which are not 
detected at all or very weak in the 21 /im sources. 

In view of the wide detection of the PAH features (and 
some flat plateau attr ibuted to HAC) an d the 30 /im feature 
in the 21 /im sources, iBuss et all (| 19931 ) suggested that the 
21 /im feature mi ght origi n ate fr om some organic molecules 
like PAHs, while iGoebel (|l993l ) argued that sulfide could 
be the carrier candidate. But the proposition of a carrier 
candidate based on its matching to these accompanying fea- 
tures is potentially problematic since they vary from source 
to source, while the profile of the 21 /im feature is rather 
universal (i.e. with the same peak wavelength, a very simi- 
lar FWHM and asymmetrical shape with a long red tail) 
indicating that the carrier of this feature should be the 
same in different sources. We also note that observationally, 
there does not appear to show any correlation between the 
strength of the 21 /xm feature with that of the 30 /im feature 
I Jiang et alj|2009h . 

For a given proposed carrier, the predicted intensity ra- 
tios of the associated features to the primary 21 /im feature 
must either be compatible with that observed (in case of 
detection) or not exceed the upper limit (in case of non- 
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detection). SiC dust with carbon impurities, once considered 
as a promising cand i date due to its abundant occurrence 
l|Speck fc Hofmeister! l2004r i. was challenged based on the 
much higher model-predicted intensity ratio of the 11.3 fim 
featu re to the 21 /im feature than observed l|jiang et al.l 
120051 ). 



4.1 Titanium-Bearing Grains 

Titanium is a relatively ra re element in the solar sys- 
tem ([Ti/H] Q ^9.77xlO- 8 , iGrevessd [l989h . and the Ti 
abundance in the meta l-poor star HP 56126 is eve n 
lower: [Ti/H]«1.3x 10" 8 (I Van Winckel fc Revniersl l2000h . 
Although some Ti-bearing grains (e.g. TiC nanoclusters and 
fullerenes coordinated with Ti atoms) have a 21 /im fea- 
ture with a close similarity to the 21 fim feature observed 
in PPNe, the deficiency of titanium would be a potential 
problem. 



4-1.1 Titanium Carbide 

IVon Helden et all |2000) argued that the 21 /im feature may 
arise from titanium carbide nanoclusters (made of 27 to 
125 atoms). The laboratory-measured spectral profiles of 
TiC nanoclusters match almost perfectly with the intrin- 
sic profile of the observed 21 fim feature, better than any 
other candidate materials, although bulk TiC dust does 
not display a strong 21 fim b and fsee Irlenning fc Mutschkel 
l200ll : IKimura fc Kaitoll2003br i. The TiC hypothesis gains its 
strength from the identification of presolar TiC grains (with 
radii ~100A) in primitive meteorites as nanometer-sized 
inclusi ons embedded in microm eter-sized presolar graphite 
grains |Bernatowicz et ahll 19961 ). However, since Ti is a rare 
element, the abundance test would be a neck to the TiC hy- 
pothesis. Indeed, three different groups have already chal- 
lenge d the TiC hypoth e sis from the Ti abu nd ance p oint of 
view ijHonv et al.l l2003l : IChigai et al.ll2003l : 12003( 1. They 
all pointed out that there may not be enough titanium to 
account for the observed strength of this feature. Here we 
add another piece of evidence against the TiC hypothesis by 
confronting it with the band strength constraint discussed in 
Sj2] Unlike previous studies, in this approach we do not need 
to know the ultraviolet/optical absorption properties of TiC. 
The experimental I R absorption spectrum o f TiC nanoclus- 
ters was measured (IVon Helden et all 20001) a nd fitted with 
a Lorentz oscillator model ( Chigai et al.ll2003h . The observa- 
tional emission spectrum of HD 56126 and the experimental 
Lorentz profile of the 21 fim band are shown in Fig.QJ,. 

The integrated mass absorption coefficient of nano TiC 
is K abs ~ 0.38 cm 3 g -1 for th e 21 fim band, as derived from 
the experimental spectrum of lVon Helden et alj ([2000 ). The 
emission temperature associated with the 21 /im carrier is 
not well constrained. 

Exposed to the stellar radiation, TiC nanocrystals, be- 
cause of their small heat capacities, will not attain an equi- 
librium temperature; i nstead, they will b e transiently heated 
by single photons (see lDraine fc Lill2001f l. The stellar atmo- 
spheric spectrum of HD 56126 peaks at A ~ 0.41 fim, i.e., a 
typical stellar photon has an energy of (hv) fa 3 eV . With 
a Debye temperature of O w 614 K (|Piersodli"996h . when 
heated by a 3 eV photon, for a TiC nano cluster as small as 



Tii4Ci3 (consisting of 3 x 3 x 3 atoms) the peak temperature 
is only T pcak w 268 K. 

We calculate from Eq. ([5} the minimum abundance re- 
quirement of Tii 4 Ci 3 dust to be [Ti/H] min w 1.97 x 10 -7 
at Td = 268 K. Since the titanium abundance in HD 56126 
is measured to be [Ti/H] t fa 1.3 x 10~ 8 (see Table [l]), the 
minimum abundance requirement (at Td = 268 K) exceeds 
the available abundance by a factor of ~ 15. 

For a given exciting photon energy (say, hv = 3eV), 
the peak temperature decreases with the nanocluster size: 
Tpcak ~ 216 K for Ti32C32 (consisting of 4 x 4 x 4 atoms) and 
Tpcak ~ 183 K for Ti72Cs3 (consisting of 5 x 5 x 5 atoms). 
Therefore, larger TiC nanoclusters would require even more 
Ti atoms to account for the observed 21 fim intensity. 

In fact, based on the available Ti abundance in 
HD 56126, nano TiC grains have to reach a temperature 
at least as high as ~ 880 K to achieve the observed emis- 
sion strength [i.e. (M d /M H ) min « (M d /M H ) maj£ ; see Fig.[2], 
while the peak temperature is only ~ 268 K even for a small 
cluster consisting of only 27 atoms. Apparently, the TiC hy- 
pothesis requires too much Ti to be viable, even under the 
most optimal condition that all Ti atoms are locked in TiC 
nanocrystals. 



4-1.2 Fullerenes Coordinated with Ti Atoms 

IKimura et ail (|2005h found that the laboratory spectra of 
large-cage carbon particles (fullerenes) coordinated with Ti 
atoms have a characteristic feature at ~ 20.3 fim closely re- 
sembling the 21 fim feature of post-AGB stars (see Fig.[lp). 
They attributed the origin of the 21 fim feature to the vibra- 
tional interaction between Ti atoms and fullerene cages. In 
order to obtain the integrated mass absorption coefficient 
re^bs on ly over t ne 21 fim band, the continuum absorption 
needs to be subtracted. Based on the continuum spectrum 
between 15-18.5 fim and 22-24 fim, we use a two order poly- 
nomial to fit the continuum underneath the 21 fim feature. 
In the following sections, unless otherwise stated, the con- 
tinuum absorption is determined in the same way. The re- 
sulting continuum-subtracted, integrated mass absorption 
coefficient for the 21 fim is K^bs ~ 0.233 cm 3 g _1 . 

Similar to TiC nanoclusters, fullerenes will also subject 
to stochastic heating in the circumstellar envelope around 
HD 56126. With a Debye temperature of 9 « 185 K0 when 
heated by a typical photon of 3 eV in HD 56126, even C60 is 
only heated to a peak temperature of T pca k ~ 89 K0 

With K^bs » 0.233 cm 3 g~ x and T pcak sa 89 K, from 
Eq. ([5j) we estimate the minimum abundance requirement 
of fullerenes with Ti atoms to be [Ti/H] min w 8.04 x 10~ 6 , 
exceeding the available Ti abundance [Ti/H] ~ 1.3 x 10 -8 
by a factor of ~618. In order not to violate the abundance 
constraint (i.e. [Ti/H] min < [Ti/H], or (Af d /Af H ) min < 
{Md/Mn) ), the dust temperature should be higher than 
400 K (see Fig.[2j| , while the peak temperature of the 
Ceo + Ti dust is only ~ 89 K. Therefore, large fullerenes co- 



2 See http://www.sesres.com/PhysicalProperties.asp 

3 For larger fullerenes, we expect a lower peak temperature T pea ^ 
and therefore require a higher minimum abundance [Ti/H] min , 
indicating a more severe Ti abundance budget shortage. 
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ordinated with Ti atoms are unlikely the carrier of the 21 /mi 
feature seen in PPNe. 



4.2 Sulfur-Bearing Grains: Silicon Disulfide 

The formation of sulfide is very likely to occur in carbon- 
rich circumstellar environ ments. As early as thirty years 
ago, lLattimer et alj (ll978T ) predicted the possible presence 
of various sulfur-bearing materials in carbon-rich systems. 
For several kinds of predicted sulfides, including MgS, FeS, 
SiS2 and CaS, the absorptio n spectra and coe fficients have 
been measured in laboratory (|Nuth et al.lll985T l. The promi- 
nent 30 ;Um dust feature seen in AGB, post-AGB and 
PNe, amounting to more than 20% of the total IR flux, 
is att ributed to MgS (|Goebel fe Moselevl Il985l ; iHonv et all 
2002). Since most of the 21 /im sources also exhibit a strong 
feature at 30 /im, it is not unreasonable to postulate that 
some sulfur-bearing grains may also contribute to the 21 /jm 
feature. Indeed, the laboratory spectrum of SiS2 displays 
a prom inent feature at ~ 22 /mi (|Nuth et al.lll985T ). IGoebel 
1 19931 ) further suggested SiS2 solids as the material respon- 
sible for the 21 /im feature. 

Based o n the laboratory spectra of IGoebel (|l993l ) and 
iNuth et al.1 (|l985l ). we obtain w 0.05 cm 3 g" 1 for the 
21 /jm band of SiS2 after subtracting the continuum (see 
Fig-HJ. In addition to the 21 /mi band, the laboratory ab- 
sorption spectrum of SiS2 also exhibits a secondary fea- 
ture at 17 /im that is never observed in the 21 /jm sources 
l|Kraus et al.l ll997T l. In order to sufficiently suppress the 
17 /jm feature so that it remains un-noticeable, the dust 
temperature needs to be < 100 K. For T<j = 100 K we es- 
timate the minimum abundance requirement of SiS2 to be 
[S/H] in ~ 9.59 x 10 -5 , exceeding the available S abun- 
dance^S/H]^ « 4.07 x 10" 6 by a factor of ~ 24. As shown 
in Fig.[2j in order for the SiS2 model to satisfy the abun- 
dance constraint (i.e. [S/H] min < [S/H] t or (M d /M H ) min < 
(M d /Mn) m „), S i S 2 du st needs to be hotter than ~200K. 

iPosch et al] (|2004h also recognized that very low dust 
temperatures (significantly lower than 100 K) would be re- 
quired to make the secondary feature at 17 /im negligible in 
strength compared to the 21 /im S1S2 feature. Although it is 
not possible to calculate the equilibrium temperature of SiS2 
due to the lack of know l edge o f its visual and near-IR optical 
constants. IPosch et al.l (|2004T ) argued that "it is hardly con- 
ceivable that SiS2 is so transparent in the visual range as to 
remain much colder than 100 K" . In our approach we actu- 
ally do not need to know the exact temperature of SiS2: (1) if 
T < 100 K - although the 17 /mi feature will be suppressed, 
one requires too much S; (2) if T > 200 K - although there 
will not be a S budget problem, the 17 /mi feature will be 
prominent. Therefore, our approach readily ruled out SiS2- 

Moreover, it is worth noting that sulphur may not be 
completely locked in SiS2. The prominent 30 /im feature, if 
indeed arising from MgS, would consume a significant por- 
tion of the sulphur available in the circu mstellar envelopes 
aroun d the 21 /im sources. According to IZhukovska fe Gaill 
( 2008), MgS has priority in the cooling sequence of sulphur- 
bearing solid compound. 



4.3 Silicon-Bearing Grains 

With four valence electrons, Si easily reacts with other atoms 
such as C, O, Fe and Mg to form various types of silicates. 
In the O-rich dust shells of evolved stars (from AGB stars to 
PNe) , both amorphous and crystalline silicates are detected 
through their numerous emission bands. Over 4000 sources 
are detected to have the most common silicate features at 
9.7 /im and 18 /im. While in the C-rich dust shells like those 
of the 21 /im sources, relatively simple Si-bearing compounds 
are formed (e.g. SiC, Si02, and SiS2). The br oad 11.3 /im 
featu re seen in five of the sixteen 21 /mi sources l|Kwok et al.l 
Il999h is identified to arise from the Si-C stretching mode 
of SiC. Being chemically active and abundant in circum- 
stellar envelopes, several Si-bearing dust species have been 
proposed to be the carrier of the 21 /mi feature. 

Silicon is an abundant element in the uni verse (its 
solar abundance is [Si/H] Q ~ 3.55 x 10 -5 , I Gre vessel 
1989), about one order of magnitude higher than sul- 
phur. Unfortunately, the Si abundance [Si/H]^ has not 
been measured for our tester HD 56126. We take the fol- 
lowing approach to roughly estimate [Si/H]^ of HD 56126: 
(1) We compile the abundance data of all 21 /im 
sources and find that fiv e 21 /im sources - IRAS 04296 , 
IRAS 22223, IRAS 23304, J Van Winckel fe Revniersl l200Ch , 
and IRAS 05113, IRAS 22272 l|Reddv et all l2002h - have 
known Si abundance. (2) We calculate the abundance ra- 
tios of Si to Ca and of Si to S for all five sources and find 
that these ratios do not vary muchQ 

(3) By assuming that the Si to Ca and Si to S abun- 
dance ratios of HD 56125 are similar to that of the other 
five 21 /im sources, we estimate [Si/H]^ m 6.31 ppm for 
HP 56126 from {Si/Ca| s a 1.55 and [Ca/H]^ « 0.18 ppm 
l|Van Winckel fe Revniersl l200Ch . or [Si/H^ « 10.7ppm 
from {Si/S} « 0.42 and [S/H]^ w 4.07 ppm. (4) We finally 
take the Si abundance of HD 56126 to be the average of that 
estimated from Ca and S, [Si/H] « 8.5 ppm. 



4.3.1 Doped- SiC 

Silicon carbide (SiC) grains with impurities were suggested 
to be the carrier of the 21 /im feature, based on labora- 
tory data that doped SiC grains e xhibit a resonance at 
~ 21 /im l|Speck fe Hofmeisterl liool ). This proposal gains 
strength from the fact that SiC is a common dust species 
in carbon-rich circumstellar envelopes. The p resenc e of cir - 
cumstellar SiC grains was first revealed by iGilral (| 19721 ). 
Nowadays, SiC is believed to be the contributor of the well- 
known IR dust feature at 11.3 /im observed in more than 



4 This may not be unexpected since S, Si, and Ca all are a ele- 
ments. Let {X/H} = log 10 [X/H] + 12 and {X/Y} = {X/H} - 
{Y/H}. For IRAS 04296, IRAS 22223, IRAS 23304, IRAS 05113, 
and IRAS 22272, the Si to Ca abundance ratios are respectively 
{Si/Ca} « 1.80, 1.57, 1.69, 1.45, and 1.49, with a mean ratio 
of ({Si/Ca}} 1.55 (and a standard deviation of a 0.12). 
The Si to S abundance ratios vary a bit more compared to 
{Si/Ca}; they are respectively {Si/S} ~ 0.70, 0.59, 0.57, 0.22, and 
0.28 for IRAS 04296, IRAS 22223, IRAS 23304, IRAS 05113, and 
IRAS 22272, with a mean ratio of ({Si/S}} RJ 0.42 (and a standard 
deviation of a K, 0.18). 
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700 carbon stars (|Kwok et al.l ll997T ). Since five of the six- 
teen 21 (im sources (i.e. about 40%) exhibit the 11.3 /im fea- 
ture, it is not unreasonable to consider both features com- 
ing from the same carrier. In particular, the distributions 
of the dust emitting at 21 /im and 11.3 /im are co-spatial 
in two 21 Mm sou rces (HD 56126 and IRAS Z02229+6208; 
iKwok et all 12002] ), which also supports the SiC hypothesis 
as the 21 /im carrier. 

However, the 21 fim feature is secondary and much 
weaker than the pri mary 11.3 /im feature in th e experimen- 
tal sp ectrum of SiC (|Speck fc Hofmeisterll20o"4h . Ijiang et al.l 
(2005) assumed a single Lorentz oscillator model for the 

21 (im feature with Q a bs(21 /im)/a treated as a free parame- 
ter (where Qabs(A) is the absorption efficiency at wavelength 
A and a is the dust size). They found that the predicted 
11.3 /xm feature is still much stronger than observed even 
with Q a bs(21 fim)/a — 10 4 cm -1 . It is unlikely to expect the 
strength of the 21 /im feature, a secondary feature of SiC 
caused by impurities, would be much stronger than that of 
the primary 1 1.3 nm feature for which <2 a b s (11.3/im)/a ~ 
1.5xl0 4 cm -4 (| Jiang et al.l2005h . Therefore, doped-SiC can- 
not produce the 21 /im feature without producing too strong 
a 11 .3 /im feature to b e consistent with observations. 

Ijiang et al.l (|2005l ) calculated the equilibrium temper- 
atures and «: a bs of doped-SiC of different grain sizes and a 
range of Q a bs(21 /im)/a values. For a-SiC, the equilibrium 
temperature is in a narrow range (~ 60-80 K), i nsensitive to 
the a ssumed Qabs(21 yxa)/a values (see Fig. 2 of ljiang et al.l 
120051 ). We see that even with Q abs (21 /im)/u = 10 4 cm" 1 
(which implies k^s m 0.25cm 3 g _1 ), the required Si abun- 
dance is at least ~8.76 x 10 -4 (for T « 80 K) and ~ 
2.31 x 10~ 3 (for T ~ 60 K), way too much compared with 
[Si/H^ « 8.5 ppm estimated for HD 56126. 

Unlike a-SiC, the equilibrium temperatures of /3-SiC 
are sensitive to Q a bs(21 /im)/u and b ecome much high er for 
smaller Q a bs(21 /xm)/a (see Fig. 3 of ljiang et al. 2005). We 
will consider cases with both large and small Q a bs(21 /im) /a. 
For Q abs (21/im)/a = 100 cm" 1 (i.e. k)£ b « 0.003 cm 3 g" 1 ), 
submicron-sized /3-SiC dust has an equilibrium temper- 
ature ~100K (see Fig. 3 of I Jiang etakl 120051 ), suggest- 
ing that one requires at least [Si/H] ; ~ 3.33 x 10~ 3 . For 
Q abs (21/im)/a = 10 4 cm- 1 (i.e. 4bs « 0.28 cm 3 g" 1 ), with 
a typical equilib rium temperature ~ 60 K (see Fig. 3 of 
I Jiang et al.ll2005h , the minimum required Si abundance is 
[Si/H] min «3.80 x 10" 3 . Therefore, for both Q-SiC and 0- 
SiC, one needs way too much Si, exceeding the Si abundance 
available in HD 56126 (~ 8.50 ppm) by a factor of two mag- 
nitudes. So it is secure to conclude that doped-SiC can not 
be the carrier of the 21 /im feature. 

4.3.2 SiC Core-SiOi Mantle Grains 

IClement et ail (|2003l ) carried out laboratory experiments 
and found that pure SiC nano-particles can be quite easily 
oxidized at their surfaces to form Si02-coated SiC grains. 
The oxidization can even reach a considerable volume frac- 
tion of the particles. The laboratory spectra of the partially 
oxidize d SiC grain s show an absorption feature at ~21- 

22 /im. iPosch et al.l l|2004 ) suggested that dust composed of 
a SiC core and a Si02 mantle may be the carrier of the 
21 /im feature. 

From the Qabs(A)/a values of IPosch et al.l (|2004l ) . we 



obtain k%£ b ss 0.21 cm 3 g _1 . The equilibrium temperatures 
of submicron-sized SiC core-Si02 mantle grai ns is ~ 104- 
126 K in the 21 /im emitting region of HD 56126 (|Posch et al.l 
120041 ). The minimum requirement of Si is ~1.00xl0 -5 - 
3.22 xlO -5 , which is slightly above the available Si abun- 
dance of [Si/H] ~ 8.50 ppm. 

However, these core-mantle grains have two prominent 
features (see Fig.[T]i) at 8.3 /im (arising from the Si02 man- 
tle) and 11.3 /im (arising from the SiC core). The 8.3 /im fea- 
ture is never seen in the 21 /im sources except three 21 /im 
sources appear to have a plateau at ~ 8 /im with an average 
FWHM of ~ 4/im which is more likely from H AC or the C-C 
stretching mode of PAHs |Kwok et al.ll20oH ). As far as the 
11.3 /im feature is concerned, we compare the flux ratio of 
the 11.3 /im band to the 21 /im band observed in HD 56126 
[F(11.3/im)/F(21/im) < 0.012] l|Honv et alj|2003h with that 
predicted from the SiC core-Si02 mantle model. It is found 
that the dust should be colder than ~ 70 K in order not to 
produce too strong a 11.3 /tm feature. However, the equi- 
librium temperature of submicron-sized SiC core-Si02 man- 
tle grains is ~ 104-126 K in the 21 /im emitting region of 
HD 56126 (|Posch et al.ll2004l ). Moreover, if the dust temper- 
ature is as low as < 70 K, 

the minimum Si abundance requirement would be 
~8.05xl0 -4 , exceeding the Si abundance available in 
HD 56126 ([Si/H]^ « 8.50 ppm) by a factor of ~ 100. 

In addition, the opacity profile of these grains has a 
shoulder in the red-wing of the 21 /im band which is not ob- 
served in the 21 /im sources. Therefore, both the secondary 
features and the 21 /im profile are discrepant with observa- 
tions, suggesting that SiC core-Si02 mantle grains are not 
a valid carrier candidate. 



4-3.3 Solid- Solution Phase of Carbon and Silicon with a 
Diamond Structure 

iKimura et ail (|2005l ) measured the IR absorption spectra of 
silicon-containing carbon films prepared by ion sputtering 
of carbon and silicon carbide mixture pellets. The carbon- 
silicon mixture film, composed of a solid-solution phase of 
carbon and silicon with a diamond structure, show signifi- 
cant absorption features a t 9.5 /im and 21 /im (see Fig.[]J). 
Thus. lKimura et alj (2005) suggested that the 21 /im feature 
observed in PPNe may arise from the solid-solution phase 
of carbon and silicon with a diamond structure. 

However, the 9.5 /im feature is far stronger than the 
21 /im feature in the carbon-silicon mixture film, but it is 
never seen in the 21 /im sources. In order to suppress the 
9.5 /im feature to such a level that the total flux emitted in 
the 9.5 /im band is less than, say 10% of that in the 21 /im 
feature, the critical dust temperature, depending on the sil- 
icon percentage, needs to be < 102 K for C-30% Si film or 
< 120 K for C-10% Si film. 

Unfortunately, there lacks a full knowledge of the op- 
tical properties of such carbon-silicon mixtures to deter- 
mine their equilibrium temperatures. We adopt the dielec- 
tric functions of diamonds which may be a reasonable ap- 
proximati on of the carbon-silicon mixtures with a diamond 
structure l|Kimura fc Kaitol2003al ) . We find the temperature 
of submicron-sized diamonds in the 21 /im-emitting region 
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of HD 56126 is ~ 148-182 Therefore, the carbon-silicon 
mixtures, if they are indeed responsible for the 21 /im fea- 
ture, would produce too strong a 9.5 /im emission feature 
to be consistent with that observed in the 21 /im feature 



4.4 Iron-Bearing Grains 

44.1 Fe 2 3 and Fe 3 4 

Iron oxide s were first suggested a s the carrier of the 21 /im 
feature bv lCoxl (|l99d ). ICox| i| 19901 ') reported the detection of 
a 21 fim emission feature in the IRAS LRS spectra of ten 
HII regions and associated it with the 21 /im band seen in 
PPNe. From the strength and the "universality" of this fea- 
ture in HII regions, iron oxides 7-Fe203 (maghemite) and 
Fe3Q4 (magnetite) were assig ned to this band. However, 
lOudmaiier fe de Winterl l)l995h later re-analyzed the IRAS 
LRS spectra of these sources and concluded that the claimed 
21 /im band was just an artifact. This was confirmed in late 
1990s when some of these HII regions were observed with 
ISO which show n o footprints of a broad 21 or 20 /im band 
l|Posch et al.ll2004 ). In fact, th e laboratory absorption spec- 
tra of iron oxides of [Co3 i|l990l ) were quite different from the 
intrinsic profile of the 21 / im feature o f PPN e which had not 
been well identified until IVolk et all l|l999h . But of course 
the spectral profile match might be improved if one consid- 
ers dust shape and clustering effects. Therefore iron oxides 
are not readily ruled out just based on the imperfect spectral 
match. 

Despite that the 21 /im feature found in HII regions 
turns out to be an artifact, iron oxides composed of cos- 
mically abundant elements and exhibiting a strong feature 
around 20 /im deserve a detailed study as being the carrier 
of the 21 /im feature seen in PPNe. In addition to 7-Fe2C>3, 
other forms of Fe2 C>3 exhibit a feature around 21 /im as well. 
iKoike et al.l (|l98ll ) indicated hematite (a-Fe203) has several 
peaks at ~ 9.2, 18, 21 and 30 /im in its absorption spectrum. 
Similarly, the absorption spectrum of Fe3C>4 shows two fea- 
ture s at ~ 17 an d 25 /im. 

ICoxl |l990l ) only presented the normalized absorption 
profiles of two iron oxides (7-Fe203 and FesCU), which 
cannot be used to calculate K a bV Based on the unpub- 
lished optical constants of Fe2C>3 and FesCU from the Jena 
grourjf], ft a b s (A) and are calculated from Mie the- 

ory and assuming a spherical radius of 0.1 /im. We obtain 
<b s w 0.09 cm 3 g _1 for Fe 2 3 and 0.07 cm 3 g -1 for Fe 3 4 . 
As shown in Fig.fJT and Fig.[TJ>', Fe2C>3 displays two strong 
bands at ~ 20.5 /im and 27.5 /im, while Fe3C>4 exhibits two 
broad features at ~ 16.5 /im and 24 /im. Apparently, FeaCU 
does not fit the 21 /im feature seen in PPNe - the model 
feature peak is neither strong nor at the right wavelength 
(see Fig.[Tg). For Fe2C>3, even with temperature fine-tuning, 



5 The Si abundance is less of an issue. From the absorption pro- 
files of iKimura et al.l d2005T) . we obtain /cjg*, « 0.32 cm 3 g -1 for 
C-50% Si film, 0.51 cm 3 g" 1 for C-30% Si film, and 0.40 cm 3 g" 1 
for C-10% Si film. If the dust temperature is respectively higher 
than ~130K, ~110K, and ~100K for C-50% Si film, C-30% 
Si film, and C-10% Si film, the minimum required Si abundance 
would not exceed that available in HD 56126. 

6 http:/ /www. astro. uni-jena.de/Laboratory/OCDB/oxsul. html#Bl 



it is very hard to suppress the 27.5 /im feature which is not 
seen in the 21/im sources. 

Moreover, Fe2C>3 and FesCU may not be able to survive 
in such a reducing environme nt like C-rich circu mstellar en- 
velopes of the 21 /im sources l|Posch et al.ll2004r ). Therefore, 
Fe203 and Fe3C>4 are unlikely responsible for the 21 /im fea- 
ture seen in PPNe. 



44.2 FeO 

iPosch et al l l|2004h pointed out that iron monoxide (FcO; 
wiistite) can survive in the C-rich reducing environment and 
proposed that nano-sized FeO dust may be the carrier of the 
21 /im feature seen in PPNe. FeO will be reduced to metal 
iron by UV photons in PNe, while it is less likely for oxygen 
to stick to iron to form FeO in AGB stars whose circum- 
stellar dust is much hotter than in PPNe. Since FeO can 
be either oxygenated to higher oxides of iron or reduced to 
iron atoms, it can survive only in a very strict physical and 
chemical environment, which is consistent with the observa- 
tional fact that the 21 /im band is rarely seen and detected 
only in PPNe that live for a short transitory period. 

We adopt the dielectr i c func tions of pure FeO measured 
by iHenning fc Mutschkei (1l997h at temperatures T = 10, 
100, 200 and 300 K to calculate the absorption spectrum of 
FeO dust. It is found that FWHM of the 21 nm band of FeO 
decreases from ~ 3.6 /im at room temperature to ~ 2.4 /im at 
T — 100 K, while the band peak A pca k shifts from ~ 19.9 /im 
to ~ 20.1 /im. At T = 100 K, both FWHM and A peak of FeO 
agree well with the observed 21 /im profile. 

Adopting the dielectric functions of FeO measured at 
T = 100 K and assuming spherical dust of radii a = 1 nm 
as sugg ested by IPosch et all (|2004T l. we calculate ^^bs ^ 
1.07 cm g for FeO (Fig.lTti shows the calculated ^ a bs 
profile which is closely reproduced by a Drude function 
with A p = 20.1/im and 7 = 2.4/im). Nano-sized FeO grains 
will be stochastic ally heated by single stellar photons (see 
iDraine fc LfeoOll). With a Debye temperature of 6 w 650 K 
l|Radwanski fc R opka 200^), the peak temperature of a 1 nm 
FeO dust is T pea k « 143 K when heated by a typical photon 
of 3 eV in HD 56126. At Td = 143 K the minimum abundance 
requirement is [Fe/H] min « 5.76 x 10 -7 . This is smaller than 
the stellar abundanc^'pe/H^ w 3.24 x 10" 6 (see Fig.0. 
Note that HD 56126 is the most iron-poor object among the 
sixteen 21 /im sources. Therefore, the FeO hypothesis will 
generally satisfy the abundance requirement. 

FeO has no noticeable secondary features in the IR, 
except a small shoulder in the blue wing of the 21 /im band 
which is not seen in the 21 /im sources. 

Stochastically-heated FeO nano dust is expected to have 
a distribution of temperatures (T < T pca k = 143 K). Al- 
though the optical properties of FeO change strongly with 
temperature, at T < 100 K they are much l ess sensitive to 
temperature (see lHenning fc M utschkc 1997). It is therefore 
reasonable to expect that FeO nano dust with a distribution 
of temperatures at T < 143 K also fits the observed 21 /im 
profile since the experimental spectrum of FeO obtained at 
T — 100 K provides the best match. 

Finally, one may ask how FeO is formed in the C-rich 
shells around the 21 /im sources where it is generally be- 
lieved that all O atoms are trapped in CO; or in other words, 
are there enough O atoms left to form FeO? Our answer is 
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"yes" . The C and O abundances (relative to H) of HD 56126 
are respectively [C/HL « 447 p pm and [O/H]^ « 468 ppm 
l|Van Winckel fe Revniera l2000l ) . A substantial fraction of 
the C atoms is required to be tied in the carbonaceous 
dust components (amorphous carbon [a-C] and HAC) to 
account for the b ulk of the IR emission of HD 56126 (see 
iHonv et al l |2003( ). Taking M a _ c « 3.6 x 10" 4 M Q and 
Mhac « 3.6 x 10" 4 M Q for the mass of amorp hous carbon 
and HAC in the dust shell around HD 56126 (|Honv et al.l 
2003), we estimate the amounts of C atoms locked in amor- 
phous carbon and HAC to be [C/H] a _ c w 68 ppm and 
[C/H] HAC fa 65 ppm, respectively^ Therefore, there will be 
at most [C/H] « 314 ppm C atoms in the gas phase to react 
with O. Assuming all gas-phase C atoms are tied with O 
atoms to form CO, there will be [O/H] « 154 ppm O atoms 
left for other O-containing molecules and dust species0Even 
assuming all Fe atoms are locked up in FeO dust, FeO only 
needs [O/H] « 3.24 ppm O atoms. 



5 SUMMARY 

We have quantitatively examined eight inorganic carrier 
candidates and one organic candidate for the 21 /im emis- 
sion band detected in sixteen PPNe, by confronting the 
abundance and accompanying IR emission features required 
or predicted by each model with that observed. We take 
HD 56126, the strongest 21 fim feature source, as a testing 
case. The principal results of this paper are the following 
(see Tabled for a summary): 

(i) Among the nine carrier candidates, five (TiC nan- 
oclusters, fullerenes coordinated with Ti atoms, SiS2, doped- 
SiC, and Si02-mantled SiC dust) are readily ruled out be- 
cause they fall short of either titanium, sulphur, or silicon. 
Even under the most optimal condition, the minimum Ti, 
S or Si abundance required to account for the observed 
21 /im feature strength still substantially exceeds the avail- 
able amount in the 21 /im sources. 

(ii) Three candidates (carbon-silicon mixtures, Fe304 and 
Fe2 0a) are ruled out because they produce strong secondary 
features (in addition to the 21 fim feature) which are not seen 
in the 21 feature sources. 

(iii) FeO nano dust, closely matching the observed 21 fim 
emission feature neither exceeds the Fe abundance budget 
nor produces undetected secondary features. There are also 

7 Here we adopt ikfjj ~ 0.44 Mg, the high end of the H mass of 
the HD 56126 shell. If we adopt the low end of Af H ~ 0.16 M®, 
the required C depletion in dust needs to be increased by a factor 
of ~ 2.8. This implies that there will be fewer C atoms in the gas 
phase to form CO and therefore more O atoms will be available 
to form FeO. We take the molecular weigh t of HAC to be /^hac ~ 
1 2.5 ( with H / (H + C) « 0.35. lHonv et al.ll2003fi . 

8 iDulevI jl980h showed that the oxidation of metallic iron is a 
very efficient process provided that enough O2 molecules are 
present. FeO i s a preferential product of low-tem perature oxi- 
dation (e.g. see lRobertslll96ll.lFehlner fe Mottlll970l, while Fe 2 3 
and Fe3 O4 are more lik e ly formed at higher temperatures (e.g. see 
iGail fe Sedrmavrlll99Sft . |Posch et al.l \2004) argued that the low- 
temperature oxidation of very small iron grains — grains composed 
of < 10 3 atoms with a size of a ^ 1 nm could result in pure FeO 
dust (instead of a large metallic iron core and a tiny FeO mantle). 



plenty of O atoms in the dust shell around HD 56126 to 
form FeO, although it is a C-rich environment. Therefore, 
nano-sized FeO seems to be a viable candidate. 

By taking an alternative approach based on abun- 
dance constraints complementary to that based on spec- 
tr al profile matching , our results are consistent with that 
of iPosch et ail fciKU: except (1) we have also ruled out two 
newly-suggested candidates (i.e. fullerenes coordinated with 
Ti atoms, carbon-silicon mixtures); (2) SiC core-Si02 man- 
tle dust is more firmly ruled out as a valid candidate; and 
(3) we prov i de fu rther support to the FeO hypothesis of 
IPosch et al.l l|2004h . Our approach does not rely on detailed 
spectral profile fitting which could be largely affected by 
dust size, shape and clustering effects. 

We have not applied the above tests to the organic can- 
didates (i.e. urea or thiourea, PAHs, and HAC) for three 
reasons: (1) the optical properties of urea or thiourea are 
not known; (2) although ten of the sixteen 21 sources 
show PAH emission features at the so-called 3.3, 6.2, 7.7, 
8.2, and 11.3 /im "unidentified IR (UIR)" bands, PAHs are 
unlikely a viable candidate since their 21 fim band strength 
is much weaker than that of the "UIR" bands; and (3) HAC 
is an ill-characterized material; its optica l properties are sen- 
sitive to the H/C and sp 2 /sp 3 ratios (see lFurton et al1ll999l ; 
I Jones et al.lfl990Tl : its 21 fim band strength (especially rel- 
ative to its vibrational bands at ~ 6-8 fim) is not well de- 
termined. To us, PAH clusters probably deserve a detailed 
investigation. 
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Figure 1. Experimental or calculated mass absorption coefficient spectra K a bsM (black solid line) in the 21 fim wavelength range of 
eight carrier candidates for the 21 (im feature seen in PPNe: (a) TiC nanoclustcrs, (b) large-cage carbon particles (fullerenes) coordinated 
with Ti atoms, (c) S1S2 dust, (d) SiC core-Si02 mantle dust, (e) carbon and silicon mixtures, (f) Fe203, (g) Fe304, and (h) FeO. Also 
shown is the astronomical emission spectrum of HD 56126 (scaled to match the K a b s peak at 21 /im; red dashed line), a proto-typical 
21 lira feature source. The blue dot-dashed line plots the continuum underneath the experimental (or calculated) 21 (im mass absorption 
spectrum. Also shown in (h) is the Drude fit (with A p = 20.1^m and 7 = 2.4/^m; green dot-dashed line) to the mass absorption profile 
K a l Ja (A) of FeO calculated from Mie theory (solid black line). The inserted panel in (h) fits the normalized, continuum-subtracted 21 ^m 
emission feature of HD 56126 (solid red line) by ft a bsM B\(T), the product of a Drude mass absorption profile (with A p = 20.1 fim and 
7 = 1.85 /jm) and a blackbody of T = 90 K (solid blue line; see Footnote-1). 
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Figure 2. f«SuL, the integrated mass absorption coefficient of the 21 (im band (with the continuum underneath the band subtracted) 
versus (■Md/MH) max = [X/H] t /nx, the maximum dust mass (relative to H) of the species containing the key element "X" (obtained 
by assuming that all X atoms are tied in that particular dust species), as denoted by cross (TiC nanoclusters) , open diamond (fullerenes 
coordinated with Ti atoms), filled circle (SiS2), open circle (doped-SiC), filled star (Si02-coated SiC), filled triangle (carbon-silicon 
mixtures), asterisk (FeO), open square (Fe203), and filled diamond (Fe304). Also plotted (as solid lines) are the minimum dust mass 
required to account for the emitted power of the 21 fim feature {M^/Mts) j = fJ-d [X/H] • /rax as a function of fc> r a range of 
temperatures. 
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Table 1. Stellar and circumstellar parameters of HD 56126 



Central Star 1 


d/kpc T cff /K 
2.4 7250 


R+/R& 
49.2 


L*/Lq 
6054 




Dust Shell 1 


Rin/cm .Rout /cm 
4.5xl0 16 9.3xl0 16 


M H /M 2 
~0. 16-0.44 






Abundance of Relevant Elements (X/H, ppm) 3 


Ti S 


Fe 


C 


O 


0.013 4.07 


3.24 


447 


168 


Total Emitted Power in 


the 21 /im Band Etot ( 


ergs -1 ) 1 


l.OxlO 36 



1 Data taken from lHonv et al . (2003). 2 Mass of the circumstellar 
shell depending on the assumed gas- to-dust ratio. 3 Data taken 
from I Van Winckel &: Revnierj (2000). 



Table 2. Tests of abundance and possible accompanying features 
for nine carrier candidates of the 21 fim feature seen in PPNe 



Candidate 


Element 


^int 

abs 


[X/H] min 


[X/H], 


Associated 


Pass (y^) or 


Material 


X 


( cm 3 g — 1 ) 


(ppm) 


(ppm) 


Features 


Fail (x) 


nano TiC 


Ti 


0.38 


0.197 (T = 268 K) 


0.013 




x 


fullerenes + Ti 


Ti 


0.23 


8.04 (T = 89 K) 


0.013 




x 


SiS 2 


S 


0.05 


95.9 (T = 100 K) 


4.07 


16.8 (im 


X 


dopcd-SiC 


Si 


0.25 


876 (T = 80 K) 


8.50 


11.3 (im 


X 


SiC + Si0 2 


Si 


0.21 


805 (T = 70 K) 


8.50 


8.3, 11.3 (Um 


X 


Si + C mixture 


Si 








9.5 fj,m 


X 


Fe 2 3 


Fe 


0.09 


3.05 (T = 175 K) 


3.24 


9.2, 18, 27.5 


X 


Fe 3 4 


Fe 


0.06 


2.92 (T = 200 K) 


3.24 


16.5,24 fim 


X 


FeO 


Fe 


1.07 


0.576 (T = 143 K) 


3.24 




V 



